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Abstract. The dissociation spectrum of the hydrogen molecular ion by short intense pulses 
of infrared light is calculated. The time-dependent Schrodinger equation is discretized and 
integrated in position and momentum space. For few-cycle pulses one can resolve vibrational 
structure that commonly arises in the experimental preparation of the molecular ion from 
the neutral molecule. We calculate the corresponding energy spectrum and analyze the 
dependence on the pulse time-delay, pulse length, and intensity of the laser for A ~ 790nm. 
We conclude that the proton spectrum is a both a sensitive probe of the vibrational dynamics 
and the laser pulse. Finally we compare our results with recent measurements of the proton 
spectrum for 55 fs pulses using a TkSapphire laser (A ~ 790nm). Integrating over the laser 
focal volume, for the intensity / ~ 3 x 10 15 W cm~ 2 , we find our results are in excellent 
agreement with these experiments. 

To be submitted to J. Phys. B 

1. Introduction 

Energy transfer mechanisms in molecules exposed to short intense laser pulses is of great 
current interest [1J. Processes such as multiple photoionization, multiphoton dissociation and 
high-order harmonic generation have been studied for a wide-range of molecular species and 
an excellent up-to-date review of the field has been provided by Posthumus 0. The hydrogen 
molecule and molecular ion represent ideal systems for a detailed understanding of ionization 
and dissociation dynamics of small molecules |3|. The simplicity of molecule means that 
it is the system of choice for theoretical studies. In general the mechanisms of dissociation 
lEfl and ionization ll20l fT6l fTvTl are very well characterised for this molecule. Nonetheless, 
the solution of three-body disintegration, incorporating spatial and temporal variations in the 
laser pulse as well as the thermal ensemble of molecular states, remains a severe challenge for 
simulation. Consequently, progress has been fairly limited in dissociative ionization spectrum 
calculations. 

Although there exists a wealth of data on the photofragment energy spectrum for the 
neutral molecule, H 2 , and its deuterated forms |2J. the data for the much simpler isolated 
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molecular ion is extremely scarce. The feasibility of experimental studies is hampered by the 
difficulty in preparing the molecular ion. In fact experiments on in intense laser fields has 
only become possible in the last 4 years ifTOl fTTl fT2l [j~3ll due to refinements in ion sources 
and charged molecular beam spectroscopy. Typically a small portion of the molecules, at a 
temperature of a few hundred Kelvin and hence predominantly in the vibrational state (v = 0), 
can be converted into the bound molecular ion by electron-impact ionization, 



The molecular ions can be extracted, cooled and collimated into a beam that can be injected 
to the focus of the laser. The photodissociation and/or photoionization processes, 



produce ions and electrons that can be collected and analysed. Owing to the small nuclear 
mass, the expansion (vibration) is very rapid - on the scale of 10 fs or less. Consequently 
the molecule relaxes extremely rapidly on the timescale of the pulse rise time. The design of 
infrared pulses that would, among other things, allow the resolution of timescale of a fraction 
of a femtosecond. 

The main aim of this paper is a direct comparison of theory and experiment for the 
proton energy spectrum produced by a short intense infrared pulse interacting with IfTOl 
fT2ll . Fragment-ion coincidence measurements that isolate the dissociation and dissociative 
ionization processes and make the comparison of theory and experiment feasible and realistic. 
The molecular dissociation dynamics are modelled using a two-state approximation and takes 
account of the laser pulse profile. Although this approach neglects the ionization channels, 
the results for dissociation energies are in excellent agreement with the observations. 

2. The model 

2.1. Physical considerations 

The characteristic vibrational and rotational time of the molecule are T vib ~ 15 fs and 
T rot ~ 170 fs, respectively. For a pulse duration 60 fs or less, the laser interaction is 
sudden on the rotational timescale. Thus one can assume a random statistical orientation 
of the molecular axis, and that the nuclei recoil along the axial direction. We assume that the 
bound ion is created by vertical transitions from the v = state of the neutral H 2 (X 1 T lg ). The 
parameters of this potential surface are, bond length _R = 1.40a.u. with vibrational constants, 
Hu>o = 0.546eV and x e = 0.0276. The v — wavepacket is promoted onto the corresponding 
molecular-ion surface if 2 l "(X 2 E 9 ) with parameters R' Q = 2a.u., D' = 2.65eV, fyjj' = 0.288eV 
and x' e = 0.0285 PTIl . The v = state is projected onto the v' manifold of states according 
to the Franck-Condon principle f34ll . The manifold of v '-states is then subjected to intense 
Ti: Sapphire light. 

In the dissociation process, momentum is conserved, for a symmetric pulse, and the 
recoiling atom and ion have equal and opposite velocity in the centre-of-mass frame. The ions 
and atoms are collected in coincidence such that dissociation and dissociative ionization can 



e" + R 2 {%, v = 0)-^e~ + e- + H+(%, v') . 



(1) 



rihv + H+( 2 S 9 , v') -> H* + H + + ri'hv 



(2) 
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be discriminated. However, depending on the method of measurement, the atomic fragments 
can arise from the entire focal volume of the laser. So a valid theoretical comparison 
must firstly establish the single molecule energy spectra, and then sum these spectra with 
the appropriate focal volume weighting for the experiment. The intensity I (p,z) of a 
focused cylindrically-symmetric laser beam is spatially Gaussian in the radial r-direction and 
Lorentzian along the axis z Q 

J 0(P> Z ) = - , ^2 eX P | 2 fl , T j ^21 [ • ( 3 ) 

where the minimum waist of the beam w and Rayleigh range Zr are given by 

w = — — , and z R = — — (4) 
irD A 

respectively. In our method, we first discuss the single-molecule response in detail, and then 

describe the orientation- averaged and focal-volume integrated results. 

The electronic coordinate, with respect to the origin at the internuclear midpoint, is 

denoted by r and the internuclear coordinate is written as R. The two lowest electronic states, 

namely the ground state ((f> g ), and A 2 Y<+ are sufficient for the study of pure 

dissociation dynamics. The corresponding adiabatic energies are E g (R) and E U (R). Then the 

two- state approximation is 

*(r, R, t) = F g (R, t)<p g (r, R) + F U (R, t)<j> v (r, R) (5) 

where F g (R,t) and F u (R,t) are the time-dependent nuclear wavef unctions. In restricting 
the electronic spectrum to the two lowest levels, the calculations will inevitably be gauge 
dependence. For low frequency fields (much less than the energy level gap) it is essential 
to use the length gauge for the laser-molecule interaction. Let us take the light as linearly 
polarized along the direction e , and denote the dipole moment between (f) u>g as n(R). Writing 
the electric field as E(t), the coupling potential is denoted by Vl(R, t) = —p(R)E(t)u, where 
u = e ■ Ft. Within the axial-recoil approximation, 9 k = cos -1 u, is the angle of ejection of the 
ion(atom) with respect to the polarization vector. Averaging over the molecular orientation is 
equivalent to averaging the projection of the electric field along the molecular axis. Since the 
dissociation rate increases rapidly with intensity over the range, I ~ 10 12 — 10 14 W cm -2 , it 
follows that the atoms and ions are ejected predominantly along the polarization direction. It 
is highly anisotropic. However at the higher intensities, / ~ 10 15 W cm -2 , the dissociation 
process begins to saturate and, as we will show, the angular distribution is broader. 

2.2. Numerical method 

Taking m = \m v to denote the reduced mass of the protons, and using atomic units, the 
two-state coupled equations are: 

i^F g (R,t) = -^^F g (R,t)+E g (R)F g (R,t)+V L (R,t)F u (R,t), (6a) 
i^- t Fu{R, t)= - ^^p F u{R> t) + E U (R)F U (R, t) + V L (R, t)F g {R, t), (6b) 
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To a very good approximation ll25ll the potential functions can be written in the form, 

E g (R) = 0.1025 [ e - 1A < R -V - 2e~ - 72 ^} , (7) 
E U {R) = 0.1025 [e- 1 - 44 ^- 2 ) + 2.22 e -°- 72 ^ 2 )] (8) 

Similarly, the dipole moment [26] can be fitted by the function ll27l 

with, p = (1 + R + R 2 /3)e~ R . For a frequency lo l and peak field strength E . This is related 
to the cycle-average intensity, / = ^c£ Eq. The time-dependence of the field can be written 
explicitly as 



E(t) = E f(t) cosco L t 



(10) 



In this paper we use the Gaussian profile that most closely models the pulses in the experiment 
of interest, 

-(41n2) ( t -^) (11) 



f(t) = exp 



T 



where f(T c ) = 1 is the maximum, and T p defines the duration of the laser pulse. 



2.3. The grids 

Discretization methods previously developed for photodissociation [32 1 and photoionization 
of molecules lfT9lfT6lfT71l can be readily applied. The discrete-variable representation (DVR) 
has proven to be a very efficient and accurate method in solving both time-independent and 
time-dependent Schrodinger equations. The Cartesian Lagrange functions [29| are given by 

/<(») = ^Vkfafakix) 



N ^ 



exp 



.2vrfc 




" 2nk 


1 N X \ 


exp 


r n x \ 



(12) 



i - (N+ l)/2 (i = 1,2,...,N), The radial coordinate R is 
discretised on this mesh Cartesian mesh between the limits -R m i n < Rj < R ma , x , with 



with the mesh points x. 
discretised 

3 e {1,2,3 



R 



N — 1, N}, such that 

(Rmax Rm 



1 N — 1 

Then the wavef unctions are expanded as: 

N 



[j ~U N+1 )]+ \( R max + Rmin) 



(13) 



F a (R,t)= J2 F AR t ,t)fi(R) 



i=i 



^ N N 

-^F CT (it^)^exp 

i=l m=l 



. 27ix m (R — Ri 
'~ N h 



[a = g,u) 



(14) 
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where h = (Rmax — Rmin)/{N — 1). Then the matrix equation corresponding to equations 
d6alandl67jl) is a dense set of 2N linear equations (er = g,u): 

N 

TijF ad (t) + V aTiij (t)F Tii (t) = iF a>j (t) (15) 

The matrix elements of the kinetic operator are given by |[29l 
/ *f (1 - i = j 

~ 1 ( 1 V-j ggg COs[7T(i-j)/JV] • / ■ (16) 

^ V AT2 s in 2 [7r(i-i)/Af] 6 ^ J- 

where the scale factor is, a = (N — l) 2 m~ 1 (R rnax — R min )~ 2 . The diagonal potentials are: 

V uu> ij(t) = 5ijE u (Rj) V gg ,ij(t) = $ijEg(Rj) (17) 

with the off-diagonal coupling: V ugyi j{t) = V gU:i j(t) = 5ijVL(Rj,t) We integrate the 
differential equations using the 18th order Arnoldi propagator as described by Peng et al 

ma. 

As is well known ll26ll . for R — > oo, the moment fi(R) — > \R, and the coupling is 
divergent in the molecular basis. This is simply a manifestation of the molecular basis and 
the use of the length gauge. It is necessary and convenient to transform to the asymptotic 
decoupled atomic eigenstates for the energy spectrum. These asymptotic states then evolve 
adiabatically at large distances and long times. However these states in turn are unbounded 
in configuration space. As shown by Keller QUI , it is possible to project the diffuse adiabatic 
asymptotic states onto a compact momentum space. More importantly, this methods produces 
the energy spectrum of the atoms that can be compared with experiment. 

Following ll28l l30l we divide the entire range of R into two regions; an internal 
region(m), where the molecular forces are significant, and an asymptotic region (as) where 
they are not. The short-range polarization potential means that the interaction region can 
be comparatively small (of the order 20 a.u.). Then the wavefunction can be partitioned as 
follows: 



where: 



F a {R,t) = F™(R,t) + F* s (R,t) (18) 

F« s (R,t) = (l-S(R))F a (R,t) (19) 
F™(R,t) = S(R)F a (R,t) (a = 1,2). 



with 



S(R) = 1 + exp 



R — i?, 



(20) 



where R s is dividing point and r s is a parameter which determines the smoothness of the 
partition. In the asymptotic region Hamiltonian is diagonalised by transforming to the 
atom+ion states xi,2- 



xT(R,t) \ = J_ ( i i \ F« s (R,t) 
X a 2 s (R,t) J V2 [ -1 1 / I F™(R,t) 



(21) 
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These states evolve into superpositions (sums over momentum) of the asymptotic states: 



1 



exp 



i(k =F A(t,t'))R - i— / dr[kTMt,r)} 2 

Jt' 



where the ion quiver momentum is 

1 * 



Ao(t,0 



dr'E{j') 



(22) 



(23) 



A symmetric or zero-area pulse is such that: A (+oo, — oo) = 0. 



2.4. Transformation to momentum space 

The key step in obtaining an energy spectrum is the projection of the numerical wavefunction 
onto the the continuum of asymptotic states (|22l) . It is efficient to do this by discretizing the 
fc-space and performing a finite Fourier transform ODI . One difficulty with such an approach 
is that the momentum shift A is a continuous function and inevitably leads to a mismatch of 
the &;-grid. In the present calculations, we use the quadrature rule for the DVR to calculate the 
Fourier transform directly. We define the shifted Fourier transform: 



Rn 



-i(k^A (t,t')R 



Xi, 2 (R,t)dR 



R 



NV2n 



exp 



N 

E 

i=l 



rn i n 

N 

E 

rn 



-i2irx m Ri/Nh 



iRmax (fc=F A (t,t')) - 



J -«[(^A (M')) -2TTX m /Nh] 
27tX r 



exp 



l(kTMt,t')) 



Nh 

2nXrr, 

Nh 



(24) 



Although this quadrature is less efficient than the fast-Fourier transform, it allows for much 
greater flexibility in choosing the k and R grids. 

At each time step, the projection onto the momentum space is made and the coherent 
momentum space wavefunction is calculated. This process is repeated and continued for 
some time after the end of the pulse. This allows the dissociating wavepacket still within the 
interaction region time to reach the asymptotic zone, and for the low energy components to 
be captured. Typically we extend the integration time to 4T C for this purpose. The A;-space 
probability density is then given by OOll : 

P k (k)= lim|xr(M) + X2 S (M)| 2 (25) 
This choice of normalization gives: 



Pfc(fc) dk 



(26) 



is the total probability for dissociation. The dissociation energy is shared equally by the 
fragments so that the proton energy is given by E p = (Am)~ l k 2 . The energy spectral density 
can then be calculated from the equation: 

P E (E P ) = (dk/dE p )P k (k) = (m p /k)P k (k). (27) 
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Figure 1. Proton energy spectrum Pe(E p ) for («' = 0) with A = 790nm and pulse 
duration T p = 55fs. The vertical lines indicate the anticipated energy of the two-photon, 
three-photon and four-photon absorption from v' = 0. 



3. Results and discussion 

3.1. Dissociation from H\ (v' = 0) at A = 790nm 

Firstly we consider photodissociation directly from v' = of H^. There are several very 
accurate proton spectra calculation for v' = dissociation by A ~ 330nm pulses for 
I ~ 10 12 — 10 14 W cm -2 |j6j|30|. We checked our calculations with these results and found 
excellent agreement in all cases. Now consider long-wavelength dissociation, we recall that 
the dissociation energy of the v' = state is D' = 2.65 eV. For A ~ 790nm, corresponding 
to the Ti:Sapphire laser, the photon energy is 1.57 eV. Therefore dissociation is necessarily a 
second or higher order process. In figure [l] we present results for the proton energy spectrum 
Pe(E p ) for a T p = 55fs pulse for a range of intensities. In these calculations, the number of 
points in _R-space iV = 512 and the number of points in momentum space N k = 2048, with 
Rmin = 0.1 a.u. and R m ax = 28.5 a.u.. For the splitting procedure we use, R s ~ 0.7R max 
and r s = 0.2. For the time-dependent propagation we use 5t = 0.01. 

For a pulse of this length, the bandwidth (FWHM) is narrow ~ 0.03 eV, and this is 
reflected in the sharply-defined proton energy (figures la, lb and lc) in the perturbation 
regime. In figure [B the relative strength of the two-, three- and four-photon dissociation 
channels change dramatically with the variation of the pulse intensity. There is a propensity 
for the high-order process. As the interaction increases there is a leftward shift of the three- 
photon peaks in frame (e) to (g). This is due to the downwards Stark shift of the ground state. 
The leading-order term of this shift is approximately linear with intensity. The broadening of 
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Figure 2. Proton energy spectrum Pe(E p ) for U^(v' = 0) with A = 790nm and 1=5.0 x 10 14 
W cm~ 2 . The figure shows the effect of increasing pulse length (a).T p = 5 fs; (b). T p = 10 
fs; (c). T p = 20 fs. The three vertical lines indicate the weak-field two-photon, three-photon 
and four-photon absorption energies. 



the peaks is evidence of the mixing of vibrational levels and shortening of the lifetime of the 
decaying state v' = 0. In the last figure (h) at the highest intensity the spectrum is dominated 
by a large broad peak around 0.7 eV where the field-free vibrational structure is destroyed. 
At intensities above 10 15 W cm -2 , the process is saturated by tunelling dissociation and the 
spectrum is dominated by low energy (less than 0.5 eV) protons [4|. This data underlines the 
difficulty in using time-independent perturbation theory to determine the ion spectrum and 
yield under these conditions. 

The second calculation (figure EJ) concentrates on the effect of pulse length for a fixed 
intensity I — 5 x 10 14 W cm -2 . For pulses shorter than 15 fs the vibrational phase is accepted 
to play a role. However it is also important to note that the bandwidth has a significant 
effect. For example, a 5 fs pulse has a bandwidth of the order 0.35 eV. This would explain 
the skewness of the peak in figure 2(a), and the presence of the 4-photon peak. However, in 
figures 2(b) and 2(c), the results resemble the narrow-band data in figure 1 . The oscillations 
in the peaks reflect wavepacket oscillations within the well and we discuss this in more detail 
later when investigating dissociation from excited vibrational states. 
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3.2. Dissociation from H 2 (v = 0) at X = 790nm 

In the ion-beam experiments, is prepared from the collisional ionization process (equation 
[T]). The H 2 (v = 0) state is projected onto a coherent superposition of the vibrational 
states, F v i, where v' = 0, 1, . . . , 17. The amplitudes C v i are given by the Franck-Condon 
factors. So that 04l in the absence of the laser the wavepackt evolves coherently 

17 

F g {R, t') = e im C vl e l ^'- lE ^' F V ,{R) (28) 

v'=0 

where the relative phases 4> v > are well defined constants, and E v > are the vibrational energies. 
The overall phase <fi(t) depends on the time of formation of the ion. According to the Franck- 
Condon factors, the largest coefficients correspond to v' = 2, 3, 4. However a small fraction of 
the population will be above the dissociation limit leading to shake-off dissociation. Dunn's 
04l estimate of this population was around 2.8%. In the present calculations, using the energy 
curves of equation ©, we obtain a value around 0.8%. In experiment the shake-off process 
is eliminated and plays no role. However, the coherence of the v' states has an important 
role. Since t' — is determined by the start of the pulse the molecules arrive in the focal 
volume with random vibrational phase. The random nature of the time delay is equivalent to 
averaging over random vibration phase. We can readily illustrate the physical consequences 
of this time delay in the proton spectrum. Suppose we choose the initial wavepacket as: 

F g (R,t' = 0) =F H2 (v = 0,t = 0) 

Then 0.8% of this wave function will naturally dissociative (shake-off) without laser 
interaction. In the calculations shown in O we present the proton energy distribution in the 
perturbation regime for a long pulse. The laser intensities are 5 x 10 9 W cm~ 2 and 5 x 10 10 W 
cm~ 2 respectively, with A = 790 nm and pulse duration T p = 55 fs. The intensity-independent 
background that peaks at zero energy is the shake-off process. The rapid oscillations in 
this background are artifacts arising from the numerical method. The lowest energy states 
emerge most slowly and require long propagation times. The major source of noise is the 
splitting procedure that creates an interference pattern in the long-wavelength wavepacket 
components. However, above 0.3 eV the interference artifacts disappear and we are confident 
in the numerical accuracy of the energy spectrum in this range. 

The 1 -photon dissociation threshold means that only levels v' > 5 appear. However 
the peak heights are modulated according to the values of both C v > and the Franck- 
Condon dissociation factors. This vibrational structure has been resolved and measured 
experimentally lfi"2l fT3l . As predicted by the perturbation theory, the probability for the 
individual peak is proportional to the laser intensity. Very recently, Serov and coworkers 
ll36l have obtained results of excellent agreement with these experimental measurements after 
careful consideration of averaging over initial ro-vibrational states and over the focus volume 
effect of the laser pulse. 

In figure © the results for high intensity show dramatic differences with perturbation 
theory. In order to compare with experiment, it is desirable to eliminate the spurious shake-off 
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Figure 3. Proton energy distribution Pk(Ek) at low laser intensities. The pulse length is 
T p = 55 fs and A = 790 nm. Full line: 1=5 x 10 9 W cm- 2 ; dashed line: 1=5 x 10 10 W cm- 2 . 
The low-energy background (dotted line) is the shake-off spectrum. The vibrational lines are 
single-photon peaks. The high intensity lines are 10 times stronger than the low intensity 
results as expected from perturbation theory. 



background. This can be achieved by projecting out the continuum or by waiting for the shake- 
off to occur and absorb the dissociation. In fact, the shake-off background is an imperceptible 
perturbation to the results at high intensity, as shown in figure |4j In practice, for intensities 
above 5 x 10 11 W cm- 2 , we have not found it necessary to remove the shake-off wavepacket 
a priori. At higher intensities the two-photon dissociation process becomes important, and in 
figure[4](b) we include vertical dotted lines to indicate the expected positions of these peaks. 
Again, at the highest intensities, the spectrum peaks move towards the lower energies although 
the spectrum also broadens. The peak heights in figure |4](c) an d (d) are similar, but the yield 
is much higher for the higher intensity. The corresponding dissociation probabilities Pd for 
each intensity are ????. As mentioned before, increasing the angle between molecule and 
polarization is equivalent to decreasing the intensity. Then figure |4| represents the angular- 
differential dissociation rates. Thus for a spatially uniform pulse of cycle-average intensity 
I — 5 X 10 14 W cm -2 , figures IH(b),(c) and (d) are respectively the yields at 9 =72, 84 and 88 
degrees. For pulses much longer than the vibrational time, the initial phase of the molecular 
wavepacket is not a significant factor. Averaging over the initial phase has little effect. In the 
next section we will consider pulses of duration comparable to the vibrational dephasing time. 
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Figure 4. Proton energy distribution at different intensities. The initial vibrational wavepacket 
is prepared from H2 (v = 0). . The pulse parameter are T p = 55 fs and A = 790 nm. The 
dashed vertical lines indicate the one-photon vibrational release with one photon absorption 
from different vibrational states and the dotted vertical line show those with two photon 
absorption. 



3.3. Vibrational phase and pulse length 

In the well-known 'pump and dump' technique, for heavy atoms within a molecule, the 
vibrational motion can be easily resolved by 100 fs pulses. However for the hydrogen atom 
the vibration is so fast (~ 15fs) that this requires pulses of 40 fs or shorter. We have already 
shown (figure El) that for H^, the vibrational phase has a role for 20fs or less. In this section 
we consider the effect on the coherent wavepacket of v ' states created from the v = state, 
and speculate on whether this can be resolved by experiments measuring the proton emission 
spectrum. We begin with a study of pulse length effects for the coherent state. In figure © the 
features of the low-intensity pulse are considered for A = 790nm and I = 5 X 10 11 W cm -2 , 
choosing the molecular vibration phases by starting the clock (pulse) at t' = 0. As noticed 
previously for v' = 0, for very short pulses the bandwidth broadening and reduced duration 
leads to a flat spectrum with low dissociation yield. Increasing the pulse duration allows 
the vibrational structure to be clearly resolved at T p = 40fs and the familiar Franck-Condon 
distribution highlighted in figure © is apparent. Also evident in figure© is the interference 
of the shake-off background with the low-energy spectrum. 
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Figure 5. Dependence of proton energy distribution on pulse length. 1=5 x 10 11 W cm -2 and 
A = 790 nm. The corresponding laser pulse length is indicated in the figure for each curve. 



At higher intensity it has been shown that the vibrational signature in the dissociation 
spectrum can be destroyed to a large extent, even for long pulses. Increasing intensity to 
I = 5 x 10 13 W cm -2 , confirms that this occurs for the coherent state (figure |6l). Firstly , the 
proton energies shift to less than leV. We note that for the 40 fs pulse, the data is very similar 
to the 55 fs pulse shown in figure[4](c). At the other end of the time scale, for 5fs pulses, the 
protons emerge with much higher energies 

That is, apart from the broadening of the peak due to bandwidth effects, the wavepacket 
favours the high frequency part of the pulse. However it is clear that in this case the pulse 
delay, or molecular phase, has an important role [37, 38 1 and ll39l l40l. In figure © we 
consider the 5 fs pulse and we present results for the wavepacket distribution in space along 
with the proton energy spectrum. The pulse is delayed by up to 20 fs with respect to the 
molecular vibration. Firstly examining the energy spectrum, there is clearly a strong variation 
with pulse delay. The 0.7eV peak that appears in figure [5] corresponding to the 5 fs peak 
with time delay fs is clearly present in figure |7J The corresponding state in R- space is 
the Gaussian wavepacket centred at R = 1.4. The short delay time allows the molecule 
to dissociate from small R values and accelerate along the S u curve. With a delay of 1 fs 
the wavepacket moves outwards and disperses, this delay severely decreases the dissociation 
yield and attenuates the dissociation energy. After 7fs delay the molecular wavepacket is 
at the outer turning point. A pulse applied at this time yields a low yield peaked around 
0.4 eV. As the delay time increases the vibrational components separate and dephase. At 
20 fs for example, the wavepacket is irregular and diffuse as the high v' components become 
evident. The dissociation from this wavepacket produces a strong low-energy signal at 0.3 eV. 
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Figure 6. Dependence of proton energy distribution on pulse length. I = 5 X 10 13 W cm 2 
and A = 790 nm. The pulse duration is indicated in the figure. 



Comparing with figure |4]for example, this energy can only be a bandwidth shifted state. That 
is, the v = 9 state dissociating through a low energy photon. In a real experiment the random 
phases of vibration mean that these details are often lost. The molecular phase (clock) can be 
synchronized by using the same laser for photoionization and dissociation 113711381 . However, 
the assumption that primary photoionization of H 2 is a vertical transition has recently been 
shown to be invalid ifTH . This is not too surprising since it is generally not a good assumption 
for infrared light at high intensity. However it is a reasonable model for UV light or electron 
impact ionization. 

Clearly the sensitivity of the proton spectrum is potentially a very useful tool for 
experimental diagnostics. The difficulty still exists in that a real laser pulse has temporal 
and spatial variations that can impede the observation of these details. For a long pulse 
T p = 40 fs, our calculations show that the kinetic energy distribution is not sensitive to the 
time-delay at the same intensity of 5 x 10 11 W cm -2 . At higher intensity of 5 x 10 13 W 
cm~ 2 , the spectrum shows only small structure effects due to pulse delay. In order to test our 
model with experiment we conclude our paper with a study in which we take into account the 
experimental parameters and integrate over the focal volume of the pulse and average over 
molecular orientation. 

i 
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Figure 7. Effect of vibration phase (time delay) for 5fs pulses. I = 5 X 10 11 W cm -2 
and A = 790 nm. (a) The i?-space probability density at the start of the pulse; (b) The 
corresponding proton energy distribution after the pulse. The different time delays are 
indicated in frame (a). 



3.4. Comparison with Experiment 

The experiment ifTOll allows the dissociative and dissociative ionization processes can be 
discriminated by coincidence time-of-flight measurements. The spectrum can be dominated 
by the large low intensity focal volume. For a theoretical model the relevant spectrum is the 
probability distribution given by 

P(I f ,E p )=J J ' 2tt pP k (I(p,z),E k )dpdz (29) 

where If is the peak intensity in the focus center and E p is the kinetic energy of the proton. 
The focus averaging over all these data with the peak intensity If = 2 X 10 15 W cm~ 2 produce 
figure [U which is in excellent agreement with the experiment for the pure dissociation peak 
ifTOl . Unfortunately, many of the structures present in the proton energy spectrum are lost after 
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Proton Kinetic Energy (eV) 

Figure 8. Proton energy distributions from theory and experiment. I = 3 x 10 15 W cm -2 
with A = 790 nm and T p = 65 fs. mo = . Full curve: Present calculation; dashed line: 
experimental measurement from the pure dissociation channel 1 10 1. 

the averaging process. The small resonant structures in the theoretical modelling correspond 
to the dissociation from different vibrational levels. At this intensity, saturation occurs at the 
centre of the focus, thus a large fraction of the dissociation fragments arise from the outer 
intensity shells. This is confirmed by our calculation in which the proton spectrum has a 
similar shape for a peak intensity 5 x 10 14 W cm~ 2 and 3 x 10 15 W cm~ 2 . Nonetheless the 
shape and ion yield of the experiment are extremely well characterised in this experiment, and 
it is highly encouraging that theory and experiment are in such good agreement. More recent 
measurement techniques, develop in the last two years, allow sections of the focal volume 
to be studied, providing and much more detailed examination of the proton spectrum. These 
new measurements will hopefully reveal some of the features of the proton energy spectrum 
described in our paper. 

4. Conclusions 

In summary, we have thoroughly investigated the dissociation dynamics of in strong laser 
fields. The vibrational- state resolved kinetic energy distributions have been calculated at 
various intensities ranging from 5 x 10 9 W cm -2 to 3 x 10 15 W cm" 2 for A = 790nm. The 
results show the sensitivity of the kinetic energy distribution on the laser intensity, molecular 
orientation, the pulse length and the pulse delay. At very high intensity, our focus-averaging 
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results reproduce the experimental measurements very accurately ifTOl . Most importantly, we 
demonstrate that molecular dissociation dynamics at rather low intensity with a short pulse 
duration would be more suitable to act as a sub-femtosecond molecular clock. 
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